Optimisation of Ionic Models to Fit Tissue Action Potentials: Application to 3D Atrial Modelling

A 3D model of atrial electrical activity has been developed with spatially heterogeneous electrophysiological properties. The atrial geometry, reconstructed from the male Visible Human dataset, included gross anatomical features such as the central and peripheral sinoatrial node (SAN), intra-atrial connections, pulmonary veins, inferior and superior vena cava, and the coronary sinus. Membrane potentials of myocytes from spontaneously active or electrically paced in vitro rabbit cardiac tissue preparations were recorded using intracellular glass microelectrodes. Action potentials of central and peripheral SAN, right and left atrial, and pulmonary vein myocytes were each fitted using a generic ionic model having three phenomenological ionic current components: one time-dependent inward, one time-dependent outward, and one leakage current. To bridge the gap between the single-cell ionic models and the gross electrical behaviour of the 3D whole-atrial model, a simplified 2D tissue disc with heterogeneous regions was optimised to arrive at parameters for each cell type under electrotonic load. Parameters were then incorporated into the 3D atrial model, which as a result exhibited a spontaneously active SAN able to rhythmically excite the atria. The tissue-based optimisation of ionic models and the modelling process outlined are generic and applicable to image-based computer reconstruction and simulation of excitable tissue

The role of infarct transmural extent in infarct extension: a computational study

Infarct extension, a process involving progressive extension of the infarct zone (IZ) into the normally perfused border zone (BZ), leads to continuous degradation of the myocardial function and adverse remodelling. Despite carrying a high risk of mortality, detailed understanding of the mechanisms leading to BZ hypoxia and infarct extension remains unexplored. In the present study, we developed a 3D truncated ellipsoidal left ventricular model incorporating realistic electromechanical properties and fibre orientation to examine the mechanical interaction among the remote, infarct and BZs in the presence of varying infarct transmural extent (TME). Localized highly abnormal systolic fibre stress was observed at the BZ, owing to the simultaneous presence of moderately increased stiffness and fibre strain at this region, caused by the mechanical tethering effect imposed by the overstretched IZ. Our simulations also demonstrated the greatest tethering effect and stress in BZ regions with fibre direction tangential to the BZ–remote zone boundary. This can be explained by the lower stiffness in the cross-fibre direction, which gave rise to a greater stretching of the IZ in this direction. The average fibre strain of the IZ, as well as the maximum stress in the sub-endocardial layer, increased steeply from 10% to 50% infarct TME, and slower thereafter. Based on our stress–strain loop analysis, we found impairment in the myocardial energy efficiency and elevated energy expenditure with increasing infarct TME, which we believe to place the BZ at further risk of hypoxia

Mediating Retinal Ganglion Cell Spike Rates Using High-Frequency Electrical Stimulation

Recent retinal studies have directed more attention to sophisticated stimulation strategies based on high-frequency (>1.0 kHz) electrical stimulation (HFS). In these studies, each retinal ganglion cell (RGC) type demonstrated a characteristic stimulus-strength-dependent response to HFS, offering the intriguing possibility of focally targeting retinal neurons to provide useful visual information by retinal prosthetics. Ionic mechanisms are known to affect the responses of electrogenic cells during electrical stimulation. However, how these mechanisms affect RGC responses is not well understood at present, particularly when applying HFS. Here, we investigate this issue via an in silico model of the RGC. We calibrate and validate the model using an in vitro retinal preparation. An RGC model based on accurate biophysics and realistic representation of cell morphology, was used to investigate how RGCs respond to HFS. The model was able to closely replicate the stimulus-strength-dependent suppression of RGC action potentials observed experimentally. Our results suggest that spike inhibition during HFS is due to local membrane hyperpolarization caused by outward membrane currents near the stimulus electrode. In addition, the extent of HFS-induced inhibition can be largely altered by the intrinsic properties of the inward sodium current. Finally, stimulus-strength-dependent suppression can be modulated by a wide range of stimulation frequencies, under generalized electrode placement conditions. In vitro experiments verified the computational modeling data. This modeling and experimental approach can be extended to further our understanding on the effects of novel stimulus strategies by simulating RGC stimulus-response profiles over a wider range of stimulation frequencies and electrode locations than have previously been explored

Revisiting QRS detection methodologies for portable, wearable, battery-operated, and wireless ECG systems

Cardiovascular diseases are the number one cause of death worldwide. Currently, portable battery-operated systems such as mobile phones with wireless ECG sensors have the potential to be used in continuous cardiac function assessment that can be easily integrated into daily life. These portable point-of-care diagnostic systems can therefore help unveil and treat cardiovascular diseases. The basis for ECG analysis is a robust detection of the prominent QRS complex, as well as other ECG signal characteristics. However, it is not clear from the literature which ECG analysis algorithms are suited for an implementation on a mobile device. We investigate current QRS detection algorithms based on three assessment criteria: 1) robustness to noise, 2) parameter choice, and 3) numerical efficiency, in order to target a universal fast-robust detector. Furthermore, existing QRS detection algorithms may provide an acceptable solution only on small segments of ECG signals, within a certain amplitude range, or amid particular types of arrhythmia and/or noise. These issues are discussed in the context of a comparison with the most conventional algorithms, followed by future recommendations for developing reliable QRS detection schemes suitable for implementation on battery-operated mobile devices.Mohamed Elgendi, Björn Eskofier, Socrates Dokos, Derek Abbot

Modelling organs, tissues, cells and devices: using Matlab and Comsol multiphysics

This book presents a theoretical and practical overview of computational modeling in bioengineering, focusing on a range of applications including electrical stimulation of neural and cardiac tissue, implantable drug delivery, cancer therapy, biomechanics, cardiovascular dynamics, as well as fluid-structure interaction for modelling of organs, tissues, cells and devices. It covers the basic principles of modeling and simulation with ordinary and partial differential equations using MATLAB and COMSOL Multiphysics numerical software. The target audience primarily comprises postgraduate students and researchers, but the book may also be beneficial for practitioners in the medical device industry

Computational model of left ventricle infarct remodelling during passive filling phase

Myocardial infarction causes part of the left ventricle muscle to not receive adequate oxygen supply and may cause muscle remodelling. Myocardium remodelling affects the blood pumping efficiency of the heart. In this paper, a model to describe myocardial remodelling during myocardial infarction is coupled with myocardium material model through multiplicative decomposition of the deformation gradient. Passive filling simulations of a normal and abnormal left ventricles are done and it was found that the muscle stiffness is increased in abnormal LV, which is indicated by 0.19% difference in the end-diastolic volume. The tissue displacement is also distributed asymmetrically in the LV, which further affects the stiffness difference. Improvements of the simulation to include FSI and full cycle simulation can further enlight the capability of the model to predict myocardium infarct remodelling

A Simplified 3D Model of Whole Heart Electrical Activity and 12-Lead ECG Generation

We present a computationally efficient three-dimensional bidomain model of torso-embedded whole heart electrical activity, with spontaneous initiation of activation in the sinoatrial node, incorporating a specialized conduction system with heterogeneous action potential morphologies throughout the heart. The simplified geometry incorporates the whole heart as a volume source, with heart cavities, lungs, and torso as passive volume conductors. We placed four surface electrodes at the limbs of the torso: , , and and six electrodes on the chest to simulate the Einthoven, Goldberger-augmented and precordial leads of a standard 12-lead system. By placing additional seven electrodes at the appropriate torso positions, we were also able to calculate the vectorcardiogram of the Frank lead system. Themodel was able to simulate realistic electrocardiogram (ECG) morphologies for the 12 standard leads, orthogonal , , and leads, as well as the vectorcardiogram under normal and pathological heart states. Thus, simplified and easy replicable 3D cardiac bidomain model offers a compromise between computational load and model complexity and can be used as an investigative tool to adjust cell, tissue, and whole heart properties, such as setting ischemic lesions or regions of myocardial infarction, to readily investigate their effects on whole ECG morphology